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Accelerated  degradation  testing  is  commonly  used  as  the  basis  to  characterize  battery  cell  performance 
over  a  range  of  stress  conditions  (e.g.,  temperatures).  Performance  is  measured  by  some  response  that 
is  assumed  to  be  related  to  the  state  of  health  of  the  cell  (e.g.,  discharge  resistance).  Often,  the  ultimate 
goal  of  such  testing  is  to  predict  cell  life  at  some  reference  stress  condition,  where  cell  life  is  defined  to 
be  the  point  in  time  where  performance  has  degraded  to  some  critical  level.  These  predictions  are  based 
on  a  degradation  model  that  expresses  the  expected  performance  level  versus  the  time  and  conditions 
under  which  a  cell  has  been  aged.  Usually,  the  degradation  model  relates  the  accumulated  degradation 
to  the  time  at  a  constant  stress  level.  The  purpose  of  this  article  is  to  present  an  alternative  framework 
for  constructing  a  degradation  model  that  focuses  on  the  degradation  rate  rather  than  the  accumulated 
degradation.  One  benefit  of  this  alternative  approach  is  that  prediction  of  cell  life  is  greatly  facilitated 
in  situations  where  the  temperature  exposure  is  not  isothermal.  This  alternative  modeling  framework  is 
illustrated  via  a  family  of  rate-based  models  and  experimental  data  acquired  during  calendar-life  testing 
of  high-power  lithium-ion  cells. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

High-power  lithium-ion  batteries  are  being  implemented  in 
various  transportation  applications  where  there  are  stringent 
weight,  energy,  power,  and  cost  requirements.  A  significant  bar¬ 
rier  that  limits  wide-spread  commercialization  of  this  technology 
is  the  requirement  of  long  battery  life,  up  to  1 5  years.  The  U.S. 
Department  of  Energy  Office  of  Vehicle  Technologies  has  sponsored 
a  variety  of  projects  to  find  solutions  to  these  and  other  barriers 
that  limit  the  commercialization  of  high-power  lithium-ion  battery 
technology.  One  related  area  of  activity  involves  the  development 
of  methods  to  accurately  predict  the  life  of  lithium-ion  batteries 
in  hybrid-electric  vehicle  (HEV)  and  plug-in  hybrid  vehicle  (PHEV) 
environments  with  a  high  level  of  precision  given  short-term  accel¬ 
erated  degradation  testing. 

Accelerated  degradation  testing  with  associated  modeling  and 
data  analysis  can  effectively  be  used  to  predict  failure-time  dis¬ 
tributions  [1].  Information  from  such  tests,  usually  obtained  at 
relatively  high  levels  of  the  accelerating  factors  (e.g.,  temper¬ 
ature),  can  be  used  to  predict  the  long-term  performance  at 
normal  use  conditions.  Typically,  accelerated  degradation  testing 
involves  exposing  cells,  in  a  collective  sense,  to  a  variety  of  thermal 


*  Corresponding  author.  Tel.:  +1  505  844  6247;  fax:  +1  505  844  9037. 
E-mail  address:  evthoma@sandia.gov  (E.V.  Thomas). 


environments.  However,  individual  cells  are  usually  exposed  to  a 
single  temperature.  During  the  time  of  exposure,  the  state  of  health 
of  each  cell  is  measured  periodically  [2  ].  The  measured  performance 
at  any  given  time  for  an  individual  cell  is  a  combination  of  effects 
that  can  be  related  to  the  technology,  to  the  unique  behavior  of  the 
individual  cell,  and  to  the  measurement  process.  Here,  we  focus 
on  an  approach  for  modeling  the  average  cell  performance  that  is 
indicative  of  a  particular  technology.  The  complete  set  of  degra¬ 
dation  data  acquired  across  all  cells  over  time  is  used  to  construct 
a  degradation  model  (e.g.,  see  [3,4])  that  expresses  the  expected 
performance  level  versus  the  time  and  temperature  under  which 
a  cell  has  been  aged.  Usually  these  models  are  cast  in  terms  of 
cumulative  degradation  at  constant  stress  conditions.  Since  cells 
are  often  deployed  in  applications  where  constant-stress  condi¬ 
tions  do  not  exist,  there  is  motivation  for  the  development  and  use 
of  rate-based  degradation  models.  By  construction,  such  models 
can  predict  cumulative  degradation  under  dynamic-stress  condi¬ 
tions.  Thus,  the  intent  of  this  article  is  to  present  an  alternative 
framework  for  modeling  degradation  that  focuses  on  the  degrada¬ 
tion  rate  rather  than  the  accumulated  degradation. 

The  remainder  of  the  paper  is  organized  as  follows.  Section 
2  discusses  cumulative  degradation  models.  Section  3  introduces 
rate-based  degradation  models  and  associated  notation.  Methods 
for  estimating  model  parameters  and  predicting  the  cumulative 
degradation  over  any  specified  temperature  profile  are  discussed. 
Section  4  illustrates  the  use  of  rate-based  models  with  existing 
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experimental  data  (obtained  under  constant  stress)  and  presents 
a  new  model  form  that  may  be  of  general  interest.  An  additional 
attribute  of  the  rate-based  framework  (not  applicable  in  the  case 
of  the  illustrative  example)  is  that  degradation  data  acquired  from 
cells  under  variable  stress  can  be  used  for  both  estimating  model 
parameters  and  model  validation.  Section  5  briefly  discusses  our 
current  thoughts  on  how  we  might  design  future  experiments  to 
take  advantage  of  the  rate-based  framework  to  use  degradation 
data  acquired  from  cells  under  variable  stress. 

2.  Cumulative  degradation  models 

The  battery  literature  contains  a  number  of  model  forms  that 
express  the  cumulative  degradation  of  the  response  of  interest  (that 
measures  state  of  health)  as  a  function  of  time  and  the  stressing  fac¬ 
tors^).  A  wide  variety  of  model  forms  have  been  proposed  (e.g.,  see 
[3]).  These  forms  can  be  mechanistic,  empirical,  or  mixed.  Aging 
temperature  is  a  common  stress  factor.  Discharge  impedance  is 
commonly  used  to  measure  state-of-health  and,  therefore,  can  be 
used  to  measure  degradation  [2].  In  such  cases,  the  expected  cumu¬ 
lative  degradation  of  a  cell  (measured  by  increasing  impedance) 
is  often  represented  by  a  model  of  the  form  /x(T ;  t)  =  a  +  k(T) 
where  /x(T ;  t)  represents  the  expected  value  of  cell  impedance  after 
exposure  to  temperature  T  for  time  t.  In  general,  we  assume  that 
/z(T ;  t)  is  an  increasing  function  of  both  t  and  T.  In  some  cases, 
we  use  the  convention  that  /x(T ;  t=0)  =  1.  This  convention  would 
pertain  to  the  case  where  the  response  variable  is  expressed  as  a 
value  that  is  relative  to  the  observed  response  at  t=  0.  For  example, 
Thomas  et  al.  [4]  considered  a  model  of  the  form, 

H(T\t)  =  1  +exp  jfo  +  jSi  •  f=}  tl/2- 

Note  that  cumulative  degradation  models  are  generally  formu¬ 
lated  in  terms  of  exposure  to  a  constant  stress  level.  Thus,  such 
models  are  often  only  used  to  predict  degradation  at,  for  exam¬ 
ple,  isothermal  conditions.  This  is  significant,  since  cells  are  often 
deployed  in  applications  where  isothermal  conditions  do  not  exist. 
While  the  instantaneous  degradation  rate  can  be  derived  by  dif¬ 
ferentiating  the  cumulative  degradation  models,  it  may  be  difficult 
to  concisely  express  the  rate  in  terms  independent  of  t.  In  such 
cases,  it  may  be  difficult  to  predict  the  level  of  degradation  in  cases 
where  the  stress  level  (temperature)  is  not  constant.  Thus,  there 
is  motivation  to  develop  rate-based  degradation  models  which, 
by  construction,  can  predict  the  cumulative  degradation  at  non- 
isothermal  conditions. 

3.  Rate-based  degradation  models 

Often,  it  can  be  more  natural  to  conceptualize  a  degradation 
model  in  terms  of  a  rate  rather  than  in  terms  of  accumulated 
degradation.  Chan  and  Meeker  [5]  describe  a  general  approach 
for  evaluating  degradation-based  reliability  that  is  based  on  a 
degradation-rate  model.  The  integration  of  the  degradation  rate 
with  respect  to  time  gives  the  cumulative  degradation  as  a  func¬ 
tion  of  the  stressing  history.  We  also  consider  models  that  can  be 
expressed  in  a  differential  form.  As  in  the  case  of  cumulative  degra¬ 
dation  models,  the  form  of  rate-based  models  can  have  a  physical, 
empirical,  or  mixed  basis.  Here,  we  restrict  rate-based  models  to 
cases  where  the  rate  of  degradation  depends  only  on  the  current 
stress  and  the  current  level  of  degradation,  but  not  the  detailed  his¬ 
tory  of  the  degradation  process.  This  is  somewhat  analogous  to  the 
“cumulative  exposure”  model  for  accelerated  life  tests  discussed  by 
Nelson  [6]  in  which  “the  model  assumes  that  the  remaining  life  of 
specimens  depends  only  on  the  current  cumulative  fraction  failed 
and  current  stress  -  regardless  how  the  fraction  accumulated.” 
Peng  and  Tseng  [7]  extends  the  cumulative  exposure  concept  to 


accelerated  degradation  testing  in  which  it  is  assumed  that  “the 
degradation  path  has  a  memoryless  property,  which  means  the  rate 
of  degradation  depends  only  on  the  current  stress,  and  not  on  the 
history  of  the  process.”  The  general  form  of  the  degradation-rate 
model  described  by  Chan  and  Meeker  [5]  also  assumes  that  the 
degradation  rate  depends  only  on  the  current  level  of  stress.  Here 
we  allow  the  rate  of  degradation  to  also  depend  also  on  the  current 
level  of  degradation.  For  example,  in  the  case  with  a  temperature  as 
the  sole  stress  factor,  d/x/dt=g(/x,  T(t)).  The  cumulative  degradation 
over  the  temperature  profile  T[0,  t]  is  therefore  given  by 


/x(T[0,  t])  = 


g(/x(r),  T(r);  oq ,  a2, . . . ,  ap)dr, 


where  g  also  involves  model  parameters  (oq,  a2,  •  •  •»  <xP)  that 
must  be  estimated.  The  model  parameters  are  generally  estimated 
using  nonlinear  regression  with  the  observed  degradation  data.  The 
observed  degradation  data  used  to  estimate  the  model  parameters 
could  be  acquired  as  a  consequence  of  non-isothermal  stress.  The 
resulting  fitted  model  can  then  be  applied  to  predict  the  degrada¬ 
tion  behavior  over  some  user-specified  temperature  profile  that  can 
be  constructed  to  match  any  particular  application  environment. 


3.  I .  Estimating  model  parameters 

The  experimental  data  are  typically  acquired  by  repeatedly 
measuring  a  number  of  cells.  To  illustrate,  consider  the  following 
notation.  Let 


Yy  =  the  jth  measurement  of  the  ith  cell  (for  i  =  1  :N  and  j  =  l:n2  ) 
ty  =  time  at  which  Yy  was  acquired,  and 

T[ 0 , tjj  ]  =  temperature  profile  experienced  by  the  zth  cell  through  ty. 

The  collection  of  data  can  be  compactly  expressed  by 

{Yu :  tu ,  T[0,  tn  ]},  {^12;  ti2,  T[ 0,  t12 f^m, ;  hn, ,  A 0,  tini  ]} 
{^21;  t2i,  T[0,  t2i]},  {Y22;  t22,  T[ 0,  t22 ]}, . . . ,  {Y2n2;  t2„2,  T[ 0,  t2n2]} 

f^ii  *  hi  ?  T[0,  hi  ]}>  { 1  h'2’  T[0,  h'2]}?  •  •  • »  {Yin. ;  tin. ,  T[0,  hnj} 

{Ym  ;  hvi ,  T[0,  hvi]},  {Vjv2 ;  t]V2,  ^[0,  tN2] },  .  •  •  ,  { YNnN‘,tNriN ,  T[0,  hvnN]} 

The  parameter  estimation  process  will  likely  involve  nonlinear 
regression.  Inherent  in  the  nonlinear  regression  process  is  a  ker¬ 
nel  function  (h)  that  predicts  Yy  by  numerically  integrating  the 
degradation  rate-based  model  over  the  temperature  profile  T[0,  ty] 
conditioned  on  the  current  set  of  model  parameter  values,  i.e., 
Yy  =  fifty ,  T[ 0,  ty  ];  ai ,  a2, . . .  oip } .  For  example,  when  using  MAT- 
LAB  (MathWorks,  Natick,  MA)  we  let  the  kernel  function  be  an 
input  to  the  nonlinear  regression  module,  nlinfit.  The  nonlinear 
regression  procedure  searches  for  the  set  of  model  parameter  val¬ 
ues  that  optimizes  some  objective  function  that  is  expressed  in 
terms  of  the  differences  between  Yy  and  Yy.  Within  the  nonlin¬ 
ear  regression  procedure,  the  kernel  function  is  computed  many 
times  for  various  sets  of  the  model  parameter  values.  Each  compu¬ 
tation  of  Yy  involves  proceeding  through  the  temperature  profile 
at  an  appropriate  time  step,  At,  using  the  recursive  relationship 
Y(/<)  =  Y(k  -  1)  +  g(Y(/<  -  1),  !(/<))  •  At  that  is  illustrated  in  Table  1. 

Estimates  of  the  model  parameters  ,  a2 , . . . ,  ap )  are  obtained 
at  the  end  of  the  iterative  nonlinear  regression  process  resulting  in 
a  fitted  rate-based  model  denoted  by 

=g(fi(t),T(t);aj,a2,...,ap). 

Note  that  the  cost  of  computing  the  kernel  function  for  rate- 
based  models  with  non-isothermal  profiles  can  be  expensive  when 


380 


E.V.  Thomas  et  al.  /  Journal  of  Power  Sources  206  (2012)  378-382 


Table  1 

Construction  of  Yy. 


8  =  t/At 

T(S ) 

Y(5) 

dm 

dt 

0 

n  o) 

Y(0) 

g(Y( 0),  T(0)) 

1 

T(  1) 

Y(0)+  .  At 

g(ni),ni)) 

2 

7(2) 

y(D+  #.At 

g(Y(2),T(2)) 

m  =  tijlAt 

T(m ) 

Y(m  l)+d^-1).A  t 

g(Y(m),  T(m)) 

compared  to  computing  the  analogous  kernel  function  for  cumula¬ 
tive  rate-based  models  with  isothermal  profiles  (where  Yy  can  be 
computed  directly). 

3.2.  Predicting  the  degradation  behavior  over  an 
application-specific  temperature  profile 

Given  the  fitted  rate-based  model,  one  can  predict  the  cumula¬ 
tive  degradation  over  an  application-specific  temperature  profile 
(ASP),  T[0,  t]ASp,  using  a  numerical  approximation  to  evaluate 
y(T[0,  t]ASP)  =  Jo  g(y(T[0,  t]„Sp),  T(r);  ai ,  a2, . . . ,  &p)dr  (e.g.,  see 
Table  1 ). 

4.  Illustrative  example 

Here  we  illustrate  the  use  of  rate-based  models  by  using  exper¬ 
imental  data  discussed  in  [4]. 

As  described  in  [4],  the  experimental  data  were  acquired  from 
high-power  SAFT  VL7P  lithium-ion  cylindrical  cells.  The  purpose 
of  the  experiment  was  to  estimate  the  mean  lifetime  of  this  cell 
technology  via  their  observed  degradation  in  performance.  Briefly, 
nine  cells  were  exposed  to  each  of  the  following  three  tempera¬ 
tures:  40  °C,  47.5  °C,  and  55  °C.  Three  additional  cells  were  exposed 
to  30  °C.  Prior  to  initiating  the  experiment  and  at  an  interval  of 
every  31.5  days  following  initiation,  the  cells  were  subjected  to 
a  reference  performance  test  (RPT)  at  30  °C  to  assess  performance 
degradation.  As  part  of  each  RPT,  the  discharge  resistance  was  mea¬ 
sured  at  a  state-of-charge  of  62%  SOC.  The  test  continued  for  about 
221  days  resulting  in  eight  RPT’s  (including  the  initial  RPT). 

Thomas  et  al.  [4]  represented  the  experimental  data  with  the 
cumulative  degradation  model, 

/4T;t)  =  l+exp{£o  +  £i-l}-tl/2’ 

where  /x(T ;  t)  represents  the  expected  relative  resistance  (resis¬ 
tance  at  time  t  divided  by  the  initial  resistance).  The  model  matched 
the  experimental  data  reasonably  well.  However,  based  on  a  boot¬ 
strap  procedure  [4],  there  was  some  indication  of  lack-of-fit.  Based 
on  this  cumulative  degradation  model,  the  expected  relative  resis¬ 
tance  (/x)  is  the  solution  to  the  equation,  (1/2 )/x2  -  /x  +  (1/2)  = 
(1/2)  exp{2  •  (fi0  +  fix  •  (1/T(r)))}dr.  For  this  particular  model 
form,  we  are  able  to  obtain  a  relatively  simple  expression  for  cumu¬ 
lative  degradation  over  a  non-constant  stress  profile.  However, 
other  model  forms  may  not  be  amenable  to  such  a  compact  expres¬ 
sion. 

We  have  found  an  alternative  model  form  that  provides  a 
closer  fit  to  the  experimental  data.  This  rate-based  model  is  of  the 
form  cfyx/dt = k(T)  ■  /n~p  with  the  initial  condition  /x(T ;  0)  =  1 ,  where 
k(T)  =  exp  {Po  +  (y0!  /T)}/(p  + 1 )  and  p  >  0. 

The  general  form  of  the  model,  d/x/dt= /<(T)  •  /x_p,  was  motivated 
by  observing  a  decreasing  rate  of  degradation  with  time  where  at  a 
particular  degradation  state  (/x),  the  rate  depends  on  temperature. 
The  Arrhenius  form  of  the  effect  of  temperature  was  assumed.  This 
model  form  is  designed  to  provide  a  flexible  family  of  rate  models 


that  depend  on  the  current  state  and  exposure  temperature  such 
that  the  rate  of  resistance  growth  decreases  with  increasing  resis¬ 
tance.  This  model,  which  might  have  application  elsewhere,  could 
apply  to  a  linear  decay  rate  where  a  reaction  was  independent  of 
the  thickness  of  a  film,  or  a  square  root  of  time  dependence  that 
mimics  resistance  of  a  film  growing  via  diffusion  limitations,  or 
something  else.  The  Arrhenius  term,  /<(T),  also  gives  this  model  form 
some  physical/chemical  basis. 

Estimates  of  the  model  parameters  could  be  obtained  using  the 
method  described  in  Section  3.1. 

However,  since  the  individual  cells  were  exposed  to  isother¬ 
mal  conditions,  one  can  use  a  simpler  approach  in  this 
case.  First,  integrate  the  rate-based  model  assuming  that  T  is 
constant.  This  yields  the  cumulative  degradation  model  /x(T, 
t)  =  [exp{fi0  +  {fixlT)}  •  t  + 1  The  model  parameters  can  be 

estimated  by  using  nonlinear  regression  with  the  data 

{Yy;ty,7)},i=l  :N,j  =  l  :  Nt. 

Following  [3],  robust  nonlinear  regression  was  used  with  the 
40  °C,  47.5  °C,  and  55  °C  data  in  order  to  estimate  the  model  param¬ 
eters.  The  estimated  model  parameters  (with  standard  errors)  are: 

fio  =  40.73  (2.2) 

fix  =-12, 194  (680),  and 

p=  11.03  (0.48). 

Confidence  intervals  for  the  model  parameters  can  be  obtained 
by  adding  and  subtracting  multiples  of  the  standard  errors  to  the 
respective  parameter  estimates.  An  approximate  95%  confidence 
interval  for  a  given  parameter  is  the  estimate  plus  or  minus  twice 
the  respective  standard  error.  For  example,  an  approximate  95% 
confidence  interval  for  fi0  is  [36.3,  45.1].  The  fitted  degradation 
model  (with  data)  is  illustrated  in  Fig.  1.  Confidence  intervals  (at 
the  95%  level)  representing  the  uncertainty  in  the  average  perfor¬ 
mance  (across  the  cell  population)  over  time  and  temperature  are 
also  provided.  One  can  compare  this  model  with  the  model  that 
was  presented  in  [4].  As  seen  in  Fig.  1,  the  current  model  gener¬ 
ally  seems  to  pass  closer  to  the  middle  of  the  observations  in  each 
time/temperature  group  than  the  model  presented  in  [4]  (compare 
with  Fig.  2  in  [4]).  While  the  rate-based  model  seems  to  under¬ 
predict  the  resistance  increase  at  30  °C,  the  apparent  discrepancy  is 
well  within  what  might  be  expected  given  a  relative  measurement 
error  of  about  1%  (see  below). 
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Date 

Fig.  2.  Hourly  temperature  profile:  Phoenix,  AZ,  Calendar  year  2010. 

Differences  between  the  observed  data  and  model  predictions 
were  assessed  in  the  context  of  the  error  model  described  in 
[4].  Briefly,  the  error  model  is  of  the  form  Var^Y^X ;  t))  ^  or|  • 
(/x(X;  t)  -  l)2  +  2  •  ol ,  where  cr|  is  the  measurement  error  and  crj 
is  a  random  cell-specific  effect.  The  estimates  of  the  error  model 
parameters  are  aj  =  3.6  x  10-3  and  a*  =  1.2  x  10-4,  respectively. 
Thus,  the  cell-specific  proportional  effect  is  estimated  to  have  a 
standard  deviation  of  about  0.06  and  the  measurement  error  is  esti¬ 
mated  to  have  a  relative  standard  deviation  of  about  0.01 .  When  the 
using  the  bootstrap  procedure  [3]  (which  incorporates  the  error 
model),  there  was  no  indication  of  lack  of  fit. 

For  this  illustrative  example,  the  assumption  of  “lack  of  mem¬ 
ory”  in  the  degradation  process  cannot  be  evaluated  due  to  the  fact 
that  cells  were  exposed  only  to  isothermal  stress.  Nevertheless, 
using  the  “lack  of  memory”  assumption,  it  is  interesting  to  apply 
the  fitted  rate  model  to  predict  the  degradation  for  a  hypotheti¬ 
cal  temperature  profile.  To  illustrate,  consider  the  hourly  ambient 
temperature  profile  of  Phoenix,  AZ  during  2010  that  is  exhibited 
in  Fig.  2  [8].  Using  the  method  discussed  in  Section  3.2  and  based 
on  the  fitted  degradation  model,  if  we  expose  cells  to  this  annual 
profile  five  consecutive  times  we  predict  the  degradation  profile 
illustrated  in  Fig.  3.  The  predicted  degradation  profile  based  on 
the  cumulative  degradation  model  from  [4]  is  also  displayed  for 


Fig.  3.  Predicted  degradation  profiles  for  rate-based  and  cumulative  degradation 
models. 


purposes  of  comparison.  Note  that  the  plateaus  in  the  predicted 
relative  resistance  relate  to  the  low-temperature  portion  of  each 
annual  cycle.  While  the  resistance  profiles  produced  by  the  two 
models  are  grossly  similar,  the  rate-based  model  begins  with  a 
gentler  increase  but  reaches  higher  levels  of  resistance.  The  unavail¬ 
ability  of  variable  stress  data  makes  it  impossible  to  rigorously 
compare  the  two  models.  Note  that  the  predicted  lifetime  of  the 
cell  would  be  determined  by  the  predicted  time  required  to  reach 
a  certain  level  of  cumulative  degradation. 

5.  Experimental  design  and  model  validation 

To  date,  due  to  the  focus  on  accumulated  degradation  models, 
most  calendar  life  experiments  involve  constant  stress.  Consider¬ 
ation  of  rate-based  models  can  motivate  a  number  of  alternative 
experimental  designs  to  assess  calendar-life.  While  one  might  con¬ 
tinue  to  expose  some  cells  to  isothermal  stress,  other  cells  could  be 
targeted  to  experience  non-isothermal  exposure.  Due  to  the  nature 
of  the  parameter  estimation  process  described  in  Section  3.1,  the 
data  acquired  from  all  cells  (independent  of  the  nature  of  temper¬ 
ature  exposure  -  isothermal  versus  non-isothermal)  could  be  used 
either  to  estimate  the  model  parameters  or  to  validate  a  model. 
It  would  be  desirable  to  subject  some  of  the  cells  to  increasing 
stress  levels  and  others  to  decreasing  stress  levels.  The  temperature 
profiles  could  be  discrete  (e.g.,  step  changes)  or  continuous  (e.g., 
sinusoidal).  One  could  also  vary  the  frequency  and  amplitude  of 
the  temperature  variations.  We  are  currently  investigating  exper¬ 
imental  designs  containing  various  combinations  of  stress  profiles 
for  studying  other  cell  technologies. 

One  could  opt  to  develop  degradation  models  with  data  acquired 
from  cells  exposed  to  isothermal  environments  and  validate  the 
models  with  data  acquired  from  cells  exposed  to  non-isothermal 
environments.  Flowever,  true  validation  of  degradation  models 
requires  that  cells  be  tested  at  stress  profiles  that  are  similar  to 
use  conditions. 

In  cases  where  the  validation  data  are  inconsistent  with  pre¬ 
dictions  from  the  developed  model,  one  needs  to  modify  the  form 
and/or  assumptions  of  the  model.  If  the  current  degradation  rate 
does  depend  on  the  detailed  past  history  of  the  degradation  process, 
then  the  degradation  measure  does  not  adequately  represent  the 
state  of  health.  Other  performance  measures  (possibly  in  combina¬ 
tion)  may  be  needed  to  represent  the  state  of  health.  For  example, 
the  growth  of  a  solid  electrolyte  interphase  (SEI)  layer  may  be 
accompanied  by  an  increase  in  resistance.  Suppose  the  resistance 
of  the  existing  SEI  layer  is  a  function  of  the  circumstances  under 
which  it  was  formed,  and,  in  turn,  affects  the  nature  of  the  layer  and 
kinetics  of  further  layer  growth.  In  such  a  case,  the  current  degra¬ 
dation  rate  would  likely  depend  on  the  detailed  prior  history  of  the 
degradation  process.  Full  characterization  of  the  state  of  health  of 
the  cell  might  require  higher-dimensional  data  (e.g.,  electrochemi¬ 
cal  impedance  spectroscopy  [9,10]).  In  such  a  case,  the  degradation 
model  would  be  more  complex  and  likely  involve  a  set  of  equations 
of  the  form: 

{lF  =Si(Mi>M2,  •••,/*<? 

=  g<j(Ml >  M 2,  •  •  • ,  Mq:  nt))|  • 

where  q  is  the  dimension  of  the  state-of-health  measure.  We  plan  to 
incorporate  such  measurements  in  our  future  degradation  studies. 

6.  Conclusion 

This  paper  provides  methodology  for  constructing  rate-based 
degradation  models.  By  construction,  such  models  can  provide 
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predictions  of  degradation  accumulated  over  variable  stress  con¬ 
ditions.  Thus,  this  methodology  can  be  beneficial  when  the  stress 
conditions  in  the  application  of  interest  are  variable  such  as  in  the 
case  of  lithium-ion  technology  used  in  transportation  applications. 
While  the  methodology  was  illustrated  with  lithium-ion  cells  under 
thermal  stress,  it  can  be  applied  more  generally  to  situations  involv¬ 
ing  other  technologies  with  other  degradation  measures  and  stress 
factors. 
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